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Abstract. We present a new *^CO(J=l-0)-line survey of the Andromeda galaxy, M31, with the highest resolution to date 
(23", or 85 pc along the major axis), observed On-lhe-Fly with the IRAM 30-m telescope. We mapped an area of about 2° x 0i5 
which was tightly sampled on a grid of 9" with a velocity resolution of 2.6kms“*. The r.m.s. noise in the velocity-integrated 
map is around 0.35 Kkms^’ on the Tnib-scale. 

Emission from the '^CO(l-O) line is detected from galactocentric radius R = 3kpc to R = 16kpc, but peaks in intensity at 
R~ 10 kpc. Some clouds are visible beyond R = 16 kpc, the farthest of them at R = 19.4kpc. 

The molecular gas traced by the (1-0) line is concentrated in narrow arm-like filaments, which often coincide with the dark 
dust lanes visible at optical wavelengths. The H i arms are broader and smoother than the molecular arms. Between R = 4 kpc 
and R = 12 kpc the brightest CO filaments and the darkest dust lanes define a two-armed spiral pattern that is well described 
by two logarithmic spirals with a constant pitch angle of 7°-8°. Except for some bridge-like structures between the arms, the 
interarm regions and the central bulge are free of emission at our sensitivity. The arm-interarm brightness ratio averaged over 
a length of 15 kpc along the western arms reaches about 20 compared to 4 for Hi at an angular resolution of 45". 

In several selected regions we also observed the '^CO(2-l)-line on a finer grid. Towards the bright CO emission in our survey 
we find normal ratios of the (2-l)-to-(l-0) line intensities which are consistent with optically thick lines and thermal excitation 
of CO. 

We compare the (velocity-integrated) intensity distribution of CO with those of Hi, FIR at 175 pm and radio continuum, and 
interpret the CO data in terms of molecular gas column densities. For a constant conversion factor Xqo, the molecular fraction 
of the neutral gas is enhanced in the spiral arms and decreases radially from 0.6 on the inner arms to 0.3 on the arms at 
R 10 kpc. We also compare the distributions of H i, H 2 and total gas with that of the cold (16 K) dust traced at /H75 pm. The 
ratios A(H l)//i 75 and (A(H l) - 1 - 2A(H2))//i75 increase by a factor of ~ 20 between the centre and R =: 14 kpc, whereas the ratio 
2A(H2 )//i 75 only increases by a factor of 4. For a constant value of Xqo, this means that either the atomic and total gas-to-dust 
ratios increase by a factor of ~ 20 or that the dust becomes colder towards larger radii. A strong variation of Xqo with radius 
seems unlikely. The observed gradients affect the cross-correlations between gas and dust. In the radial range R = 8-14 kpc 
total gas and cold dust are well correlated; molecular gas is better correlated with cold dust than atomic gas. At smaller radii no 
significant correlations between gas and dust are found. 

The mass of the molecular gas in M 31 within a radius of 18 kpc is Af(H 2 ) = 3.6 x 10* Mq at the adopted distance of 780 kpc. 

This is 12% of the total neutral gas mass within this radius and 7% of the total neutral gas mass in M 31. 
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1. Introduction 

Star formation and spiral structure in galaxies require the cou¬ 
pling of processes operating on linear scales so different that 
they are hard to study in a single galaxy. The small structures 
are difficult to observe in external galaxies, whereas large struc¬ 
tures, due to distance ambiguities, are hard to see in the Milky 
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Way. Single-dish telescopes were used to survey CO in galaxies 
(e.g. Nakano et al. nWTl Braine et al. IT^ Young et al. 11995} 
but with limited angular resolution. Molecular spiral arms were 
barely resolved in these surveys even in the nearest galaxies 
(e.g. Koper et al. 119911 Garcia-Burillo et al. 119931 Loinard et 
al. ll996l Heyer et al. l2004l l. Only in the Magellanic Clouds the 
single-dish surveys resolved giant molecular clouds (Israel et 
al.[T993}. More recently mm-wave interferometer surveys like 
the BIMA SONG (Regan et al. 1200 U gave vastly improved 
data on nearby galaxies like M51, resolving molecular arms 
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into cloud complexes. This instrument was also used for an all¬ 
disk survey of M 33, about ten times closer to us than M 51, in 
which individual molecular clouds are recognized (Engargiola 
et al. l2003> . The IRAM Plateau de Bure interferometer has re¬ 
solved molecular clouds in M 31 into components (Neininger 
et al. l2()()()at enabling close comparisons with molecular clouds 
in the Milky Way. 

The nearest large spiral is the Andromeda Nebula, M 31. Its 
distance of D = 0.78 + 0.04 Mpc (Stanek & Garnavich ll998t 
ranks among the best known for any galactic or extragalactic 
nebula. The accuracy of this distance allows us to derive ac¬ 
curate luminosities and masses. At this distance 1' along the 
major axis corresponds to 227 + 12 pc. The large inclination of 
M31, i = 77°5, degrades the resolution along the minor axis 
by a factor of 4.6, but has the advantage of yielding accurate 
in-plane velocities. M 3rs proximity gives us the chance to see 
many details of the distribution and kinematics of the gas as 
well as the relation of the gas to the spiral structure and to star 
formation. 

The contents of stars, dust, and atomic gas in M 31 are well 
known. The whole galaxy has been mapped in the 21 cm line 
of Hi with 24" x 36" resolution by Brinks & Shane (119841 
hereafter B&S) and its northeastern half with 10" resolution by 
Braun (I1990t . It has been entirely mapped in the mid and far 
infrared by the IRAS, ISO and Spitzer satellites (see Haas et 
al. ( fT^ and Schmidtobreick et al. ( EUnOt for the ISOPHOT 
map at 175 fim, and Gordon et al. ll2004t for the MIPS maps 
at 20 jum, 60and 160 fim). Furthermore, M 31 was partially 
mapped with ISOCAM (5.1-16.5 fim) at 6" resolution (see e.g. 
Pagani et al. 119991 . Comparisons of the emission in different 
wavelength ranges - like UV, optical. Hi, FIR (160/im and 
175 ijm) and radio continuum emissions - have also been re¬ 
ported (Foinard et al. 119991 Pagani et al. 119991 Keel 120001 
Feaueux 120001 Nieten et al. 120001 Berkhuijsen et al. 120001 
Gordon et al. 120041 . 

So far, the situation was not as favourable as to the molecu¬ 
lar gas. Prior to ours, the only complete CO survey of M 31 was 
made with a 1.2-m diameter telescope and had a resolution of 
8!7 (Koper et al. 119911 Dame et al. 119931 . More recently, a sur¬ 
vey of the southwestern half, made at an angular resolution of 
1' with the FCRAO 14-m telescope, was published by Foinard 
et al. (1199611199^ . The latter authors (119991 their Table 2) give 
a nearly complete overview of previous CO observations of 
M 31. Foinard et al. ( I1999> and Heyer et al. (12()()()t found many 
similarities, but also clear differences, between properties of 
the molecular gas in M 31 and those in the Milky Way. 

Our survey, made with the IRAM 30-m telescope in the 
'^CO(J=1-0) line, has a resolution of 23" corresponding to 
85 pc along the major axis. It is much more sensitive than 
the previous surveys and detects all clouds with J Tmbdv ^ 
1 Kkms * (= 3x rms noise). In this article we present the CO 
distribution in the bright disk of the galaxy. We derive some im¬ 
portant basic results using simple assumptions, e.g., a constant 
conversion factor Aco from CO intensity to molecular column 
density. We discuss the spiral-arm structure of the neutral gas 
and the arm-interam brightness contrast in Sect. 3. In addition 
to the *^CO(J=1-0) line, we have observed several selected 


areas covering bright arm segments, in the ^^CO(J=2-l) line 
with high sensitivity; we discuss the line ratios in Sect. 3.3. In 
a previous publication, based on one third of the present data 
(Neininger et al. ll998> . we reported a tight correlation between 
the CO sources and the dark dust lanes. In Sect. 4 we return 
to this point and compare the CO distribution with those of 
Hi, fir (175/im) and A20 cm radio continuum. Radial pro- 
hles of the various constituents are discussed in Sect. 4.1 and 
correlations between CO, Hi and FIR (175/im) in Sect. 4.2. 
In Sect. 4.3 we derive the molecular and total gas mass. The 
CO velocity held is described in Sect. 5. Our results are sum¬ 
marized in Sect. 6. Preliminary reports on this survey were 
given by Guelin et al. (|2000j, Neininger et al. (|T99i|2000B 
and Nieten et al. (|2000j. 

2. Observations 

Our survey was carried out with the IRAM 30-m telescope 
between November 1995 and August 2001. The observations 
were made On-the-Fly in two steps; in a hrst step a held typ¬ 
ically 18' X 18' in size was scanned back and forth in the di¬ 
rection parallel to M3rs minor axis, Y, at a speed of 4" s“^ 
The successive scans were spaced by 9" in the orthogonal X 
direction. At the beginning and at the end of each scan, a ref¬ 
erence position, located 30' or 45' away from the major axis' 
and free of CO or Hi emission, was observed for 30 s. Every 
1-2 hours the telescope pointing was checked on planets and 
nearby quasars. A second reference position, located within 
M 31 and showing strong CO emission, was observed for cali¬ 
bration purposes (see below). The telescope focus was checked 
several times a day, in particular after sunrise and sunset. In 
a second step the observations were repeated by scanning the 
same held in the orthogonal direction, parallel to the major 
axis. The data recorded by the backends were read every sec¬ 
ond of time, so that the data cube obtained by combining the 
two orthogonal maps was fully sampled on a 9" grid. 

The reduction procedure was described in some detail by 
Neininger et al. (1200061 . After calibration (see below), sub¬ 
traction of the off-source reference spectrum and a baseline 
for individual spectra in a map, two orthogonal maps were 
combined using “basket-weaving”, the de-striping technique of 
Emerson & Grave ( 119881 . This code was adapted to work on 
two orthogonal channel maps before averaging them (Hoernes 
1 19971 . Examples of the technique are shown in Neininger et al. 
(I2000bt . Hoernes ( 119971 and Emerson & Grave (EMHll. In this 
process, the CO map was smoothed from 21" to 23" FWHM. 

The maps shown in Fig. ^are the combination of 12 in¬ 
dividual helds, listed in Tabled Each held is fully sampled 
and all helds together contain nearly 1.7 million spectra (be¬ 
fore gridding) obtained in about 500 hours of effective observ¬ 
ing time. We present individual spectra for selected regions in 

Fig.El 


' Throughout this paper we use the M31-hxed coordinate system 
of Baade & Arp <19641 : centre coordinates RA = 0''40"'0(7.3, DEC = 
41°00'03" (1950.0) (Dressel & Condon liWSl . X along the major axis 
(PA = 3TJ), positive to the NE, Y along the minor axis, positive to 
the SE. All velocities are in the LSR frame. 


























































Ch. Nieten et al.: Molecular gas in the Andromeda galaxy 


3 



30 10 0 -10 -3& 30 LO 0 -10 -30 

T offset [nromm] Y ofbrt [flrcmiii] 

Fig. 1. (a) The velocity-integrated intensity distribution of the '^CO(l-O) spectrum, /i_o = / Tinb(C(9i_o)dv, observed with the 
IRAM 30-m telescope. The X and Y coordinates are taken along the major and minor axis, respectively; the position angle and 
centre position are taken from Dressel & Condon J 19761 - see footnote 1). The dashed line marks the border of the area surveyed 
which is about one degree squared. The colour scale starts at l.SKkms"'. Most structures visible in the map, including weak 
ones, correspond to line profiles detected with signal-to-noise ratios > 4. The mean r.m.s. noise in the integrated intensities 
is 0.35Kkms ', but the noise varies across the map (see Table [Q. (b) The velocity field (first moment) as traced by the CO 
emission. Velocities are LSR-values. 
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When observing the CO(l-O) line, the 30-m telescope al¬ 
lows simultaneous observations of two polarizations of the (2- 
1) line. The (2-1) line, however, is weaker than the (1-0) line 
and the integration time per beam is about 4x smaller because 
the beam area is 4x smaller. Moreover, the receivers are nois¬ 
ier and the sky opacity is higher at 230 GHz (2-1) than at 
115 GHz (1-0). This made it difficult to detect the (2-l)-line 
emission in our survey, except from the brightest clouds. In or¬ 
der to improve the signal-to-noise ratio and the sampling for 
(2-1) we re-observed several rectangular regions, 3'^' wide 
by 3'-12' long, with half the scanning velocity and twice the 
sampling, when the zenith opacity at 230 GHz was favourably 
low (< 0.2). We discuss some general results of these obser¬ 
vations in Sect. 3.3. In order to confirm the reliability of the 
OTF method, to integrate some emission-free positions to a 
lower noise level, and to check several apparent discrepancies 
with previous CO observations we re-observed about 200 posi¬ 
tions, located inside as well as outside the arms, in the position¬ 
switching mode with integration times of 10-30 min. Some of 
these results are used in Sect. 3.2.^ 

We used two SIS receivers with orthogonal polarizations to 
observe the (1-0) line and a similar system to observe the (2-1) 
line. The receiver temperatures in the standard reference plane 
(before the polarization splitter) were close to 90 K (SSB) at 
115 GHz at the beginning of our survey (fields 1-3, see Tablellli 
and close to 50 K at the end (fields 8-11). After addition of the 
atmospheric contribution, the system temperature was between 
200 and 400 K at both 115 GHz and at 230 GHz. The backends 
consisted of two 512 x 1 MHz filterbanks at 115 GHz and of 
2 autocorrelators with resolutions of 0.8 MHz and total band¬ 
width of 320 MHz at 230 GHz. A channel width of 1 MHz 
corresponds to a velocity resolution of 2.6 km s * for the (1-0) 
line and of 1.3 km s“' for the (2-1) line. 

The standard calibration at the 30-m telescope is equiva¬ 
lent to the chopper-wheel calibration method for observations 
at millimeter wavelengths and gives the antenna temperature 
T*^, corrected for atmospheric losses and forward efficiency, 
Ff-fi (for details see Downes ri989t . The main-beam brightness 
temperature, can then be calculated from 

Tmb = F’eff/BeffT; , (1) 

where Beff is the main-beam efficiency. Over the years of this 
project, Beff and in particular improved significantly due to 
upgrades of the 30-m telescope, the most important of which 
occurred in July 1997. We used the standard values for these 
efficiencies for the last quarter of 1997 derived by the obser¬ 
vatory staff, for Beff from antenna tippings and for Beff from 
observations of Mars and Uranus: Beff/Beff = 0.92/0.80 at 
115 GHz and 0.80/0.50 at 230 GHz. Day-to-day variations 
were monitored by observing two reference positions: X,Y - 
-30H8, 4f26 in the south and X,Y - 18f80, U90 in the north. 
Normalization of the line integrals at these positions to the 
same values ensured the same calibration parameters for the 
entire survey. The uncertainty in the final B^b scales is about 
15%. 

^ The full results of the position-switched observations will be pre¬ 
sented elsewhere. 


A corollary of the high beam efficiency of the 30-m tele¬ 
scope at 115 GHz is the low error beam. Greve et al. ifT^ 
found that the far-beam pattern can be described by 3 error 
beams with half-power widths of 5', 7' and 80', respectively. 
The first two could pick up signals from regions with similar 
radial velocities as observed in the main beam. As they con¬ 
tribute only a few percent to the main-beam signals, they hardly 
affect the observed spectra. The contribution of the third error 
beam, which is larger than the disk of M 31, is negligible. 

The effective integration time per 23" beam for the large 
'^CO(l-O) map was 64 s, yielding a r.m.s. noise of ^ 33 mK 
per 1 MHz channel (Bnib scale) in the southern fields (1-4) and 
^ 25 mK per 1 MHz in the northern ones (see Table Q. The 
corresponding values for the small maps are 21", 173 s and 
15 mK for the (1-0) line, and 12", 57 s and 35 mK per 1 MHz 
channel for the (2-1) line. 

The distribution of velocity-integrated CO-line intensities, 
7i-o = f 7’inb(CGi_o)dv, is shown in Fig. The r.m.s. noise 
in the velocity-integrated emission varies between the fields, 
but is typically about 0.38Kkms ' in the southern part and 
about 0.33 Kkm s ' in the northern half of the survey or 0.35 + 
O.lOKkms ' for the total map (see Table We note that the 
sensitivity of our survey to point-like and extended sources ex¬ 
ceeds that of Loinard et al. (tT9^ by factors > 8 and > 1.5, 
respectively. 

3. The CO brightness distribution 

The survey of the Andromeda galaxy presented in this paper is 
the largest and most detailed molecular-line survey ever made 
of an extragalactic object. Most of the emission in Fig. □ ap¬ 
pears concentrated to radii between 3 and 12 kpc and occurs 
as long and narrow filaments that strongly suggest a spiral arm 
structure. In addition, we see a number of scattered CO clouds 
of weak intensity between the spiral arms; sometimes these 
form bridges. The inner arms at radii near 5 kpc are remark¬ 
ably bright, especially in the NE half of the galaxy. 

Close to the centre of M 31 CO emission is very weak. 
Melchior et al. J2000t have found CO emission of about = 
20 mK in a dark-cloud complex located at a distance of 350 pc 
from the centre using position switching and long integrations. 
The integration time per beam was about 6 hours and the r.m.s. 
noise 2 mK per 3 km s '-wide channel, which is much better 
than the noise of the present survey. 

Only a few scattered clouds are visible at large distances 
from the centre. The most distant cloud in this survey was 
found at a deprojected radius of 19.4 kpc. The spectrum of this 
cloud is shown in Fig. Et The central velocity of the emission 
fits well to the velocity field of M 31 (see Fig. m. The cloud is 
located in the outermost part of the spiral arm A-N modelled 
by Braun J1991t near several H ii regions. 

In the following subsections we model the spiral pattern 
and analyze the arm-interarm contrast. 

3.1. The Structure of molecular spiral arms 

Although M 31 was classified from the beginning by Hubble 
il929l as a Sb type spiral, attempts to draw its spiral pattern 
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Table 1. Data on observed fields of M 31 


Field 

# 

Quadrant 

inM31 

A-extent 

(') 

T-extent 

(') 

Number 
of spectra 

Mean r.m.s. noise 
in /i_o per channel 

(Kkrns-*) (mK)‘> 

1 

S 

-48.5 ... - 

33.3 

-10.1 ... 

14.0 

162980 

0.45 

34 

2 

SW 

-34.0 ... - 

13.9 

-18.1 ... 

-2.0 

225122 

0.38 

37 

3 

SE 

-34.0...- 

13.9 

2.0... 

16.1 

115765 

0.35 

34 

4 

SW 

-14.2 ... 

5.1 

-17.0 ... 

-3.0 

128027 

0.44 

33 

5 

NE 

21.0... 

39.3 

-2.0... 

16.1 

204190 

0.46 

32 

6 

NW 

21.0 ... 

39.1 

-19.6 ... 

-1.5 

152630 

0.25 

22 

7 

NE 

39.0... 

55.9 

-2.0... 

15.1 

160045 

0.27 

21 

8 

NW 

38.9... 

56.1 

-19.6... 

-1.5 

117198 

0.29 

27 

9 

NE 

3.0... 

21.1 

-2.0... 

18.0 

164760 

0.32 

25 

10 

SE 

-14.1 ... 

3.2 

-3.0... 

17.0 

128534 

0.27 

27 

11 

NW 

4.8 ... 

21.1 

-19.0... 

-1.0 

128156 

0.33 

24 

12 

NW 

2.8 ... 

5.1 

-3.0... 

-1.8 

666 

0.36 

26 


NE 

A> 0 


T > 0 

528995 

0.362> 

262> 


NW 

A> 0 


T < 0 

398650 

0.29 

24 


SE 

A< 0 


T > 0 

339371 

0.35 

31 


SW 

A < 0 


y < 0 

421057 

0.41 

35 


Total 

3512 square arcminutes 


1688073 

0.35 

29 


Temperatures are on T^ij-scale, where Tmi, = 1.157^. Weighted means 


have been mostly inconclusive. For example, Baade (119631 
used the young stars and Hodge ( fPT^ the open star clusters as 
spiral arm tracers. The results have been summarized by Hodge 
( I1981bt . Due to absorption of optical light by dust-rich lanes, 
the presence of a bright optical bulge and the lack of H i radio 
line emission in the inner disk, there was no consensus even on 
whether the arms are trailing or leading. More recently, Braun 
( I1991> proposed a non-planar trailing two-armed spiral pattern 
with varying pitch angle on the basis of the H i interferometric 
data. His model accounts fairly well for the outer H i arms, but 
it does not trace the innermost structures. Based on their anal¬ 
ysis of the 175 /tm dust emission, Haas et al. (1199811 suggested 
that M 31 may be closer to a ring galaxy than to a spiral. 

The CO emission, which traces the dense molecular gas, is 
better suited to determine M 31 ’s spiral pattern because the CO 
arms are thinner than the Hi arms, they are less patchy than 
the Ha arms and not affected by absorption. As is discussed 
below, the CO arm-interarm contrast is also much larger than 
in H I. Furthermore, the linear resolution of the CO survey is 
sufficient to distinguish neighbouring arms even on the minor 
axis. 

In order to determine the spiral pattern in an objective way, 
we decomposed the CO map into 170 individual ‘clouds’, us¬ 
ing the CLUMPFIND analysis program (Williams et al. ll994t . 
CLUMPFIND provided the position of ‘clouds’, deprojected 
them (assuming an inclination of the molecular disk of i = 
7T.5) and transformed them into polar coordinates. The result¬ 
ing data points, plotted in semi-logarithmic coordinates, were 
then least-squares fitted with two straight lines representing 
two simple logarithmic spirals (see Figs. 1^ and b). The fits 
led to very similar pitch angles for the two spirals, namely 


i/r = 8'’ + D and 1° ±1°, and to a phase shift of roughly 200'’. 
The derived pitch angles are in good agreement with the value 
of i// - 6°7 for the mean of the 2 spiral arms observed in H i 
(Braun I1991t and with the mean value of 7°A for the optical 
spiral arms (Baade & Arn 119641 . The oscillations visible in 
Fig-Eh indicate variations in pitch angle which could be partly 
due to variations in inclination angle observed by Braun (I1991> . 

We note that these geometrical spiral fits required cen¬ 
tre positions shifted to the NE with respect to the nucleus 
of M31, i.e. Xc,Yc - 7'.7, O'A for the full-line spiral and 
Xc, Yc - 5'8, 2(3 for the dashed spiral (Figs.|^ and b). Such 
a displacement is not surprising as the bulge of M 31 is likely 
to host a bar (Stark & Binnev [19941 . which makes it difficult 
to trace the spiral arms down to the nucleus. This displacement 
and the fact that on the minor axis the outer CO arms lie slightly 
outside the fitted spirals may also indicate that the molecular 
spiral arms have inclinations that differ from the adopted in¬ 
clination of 77°.5 of the main plane of M 31, like the Hi arms 
fBraun [T991t . For a detailed comparison with the Hi arms a 
more complete analysis of the molecular arm structure should 
be made, which also includes the velocity structure of the arms. 
This could be the subject of a future study. 

For comparison, we have tried to represent the CO emission 
atR- 8-12 kpc by a circular ring, as suggested by Haas et al. 
(179^ . The result of the fit is shown in Fig. Et. Even at these 
radii the spiral pattern appears to be a better description of the 
CO distribution than a ring-like structure. 

The derived spiral pattern is a good fit to the CO emission 
up to about 12 kpc from the centre, except for the region of the 
giant stellar association NGC 206 (X, Y - -40', +2'), the ve¬ 
locity of which strongly deviates from the normal rotation ve- 
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Fig. 2. Sample of spectra extracted from the data cube of the ^^CO(l-O) survey of M31 (full lines). At some positions also 
the '^CO(2-l) spectrum is shown (dotted lines, shifted down by 0.1 K for clarity). The velocity resolution is 2.6 km s ' for 
the (1-0) line and 1.3 km s ' for the (2-1) line. Velocities are in the LSR frame, (a) Weak emission at V, T = 46(1, -15(6 
which corresponds to a deprojected distance from the centre of 19.4 kpc. (b) Spectrum with linewidth < 5 km s ' at V, T = 
-46(08, 9(35. (c) Broad line of about 40 km s ' close to the minor axis at V, T = -0(33, -3(49. The (2-1) spectrum is smoothed 
to 9.9kms '. (d) Two peaks separated by about 40kms ' at V, T = 3(02, -4(97. (e) Emission from the dust cloud D84 at 
X,Y - -16(528, -8(728; strongest emission on 10-kpc arm. (f) Emission at the position V, y = -17(8, -3(85 near the dust cloud 
D 153. (g) Emission with three velocity components near the dust cloud D 47 at V, T = -22(48, -7(53. (h) Double peak emission 
from the inner arm, X,Y — -7(57, 5(115. (i) Emission with three components near the dust cloud D 39, X,Y - -41(9, 8(54. (j) 
Emission at the position X,Y - 24(85, 1 (62 near the dust cloud D 615. 


locity. In particular, it accounts for the splitting of the CO emis¬ 
sion along the eastern minor axis into two distinct arms. These 
arms are clearly seen in the 8 /rm emission from hot dust (MSX 
satellite map, Moshir et al. 119991 ISOCAM map of the central 
bulge of Willaime et al. 120011 see Eig.0 below) and are also 
predicted by Braun’s model although they are only marginally 
visible in H i and Ha emissions. 

At radii ^ 12 kpc the fitted spiral pattern starts to devi¬ 
ate from the observed spiral. There the CO emission becomes 
fainter and the filaments more difficult to trace. We note, how¬ 
ever, that the ISO 175 /urn map of Haas et al. lfT9^ and the Ha 


image of Walterbos J2000t show a clear arm-like structure that 
seems to lie on the extension of the fitted spiral to the NE, far 
beyond our CO image. 

3.2. Arm-interarm contrast 

To derive the apparent arm width and the arm-interarm bright¬ 
ness intensity ratio, we divided M 31 ’s disk into 4 quadrants of 
approximately equal size. Eor each quadrant, we generated a 
family of spiral segments from a linear combination of the spi¬ 
rals that best fit the inner and outer CO arms. Eigure|^ shows a 
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Fig. 3. Spiral structure in M 31. (a) Two logarithmic spirals (solid and dashed lines) htted through the distribution of CO ‘clouds’ 
in the radial range R - 4-12 kpc. They are overlayed onto the CO brightness distribution (contours), (b) Least-squares hts in the 
\nR - (j) plane yielding pitch angles of 8'’ and 1° for the solid and dashed spirals in (a), respectively, (p is the azimuthal angle in 
the plane of M31 counted clockwise from the eastern minor axis, (c) Ring htted to the distribution of ‘clouds’ at R = 9-12 kpc 
(solid line) superposed onto the CO-brightness distribution (contours). 


sample of the spiral segments generated for the NW quadrant, 
superposed onto the CO map smoothed to a resolution of 45". 
We then calculated the average CO intensity along each spiral 
segment. Figurel^i shows these averages for the 4 quadrants as 
functions of the spiral segment radii at their intersection with 
the major or minor axis. In this analysis, we have discarded 
the highly perturbed region surrounding NGC 206 in the SW 
quadrant. 

The CO peak-to-dip ratio reaches a maximum > 20 in 
all quadrants, showing that the arm-interarm intensity contrast 
is high almost everywhere. The shallow dip at /? = 11 kpc 
in the SE quadrant reflects that the two brightest spiral arms 
{R -9 kpc and 12 kpc) are only partly resolved in this quadrant 
(see Fig.[0. The shoulders visible on Fig. HatR = 9 kpc (NE 
quadrant) and at 10 kpc (SW quadrant) come from two weak, 
nearly symmetrical arm-like substructures that are clearly vis¬ 
ible on Fig.n(e-g from X,Y - 27', 8' to X,Y = 40', 2' in the 
NE quadrant). 

That CO emission is very low between the arms, outside 
these substructures, is supported by the more than one hundred 
pointed observations we did in the interarm regions, in particu¬ 
lar in the direction of several weak and extended interarm fea¬ 
tures that appear on the map of Loinard et al. (I1999t . Our obser¬ 


vations failed to reveal any emission down to levels 3-5 times 
lower than the sensitivity of our OTF survey. Similarly, a strip 
of 8' length centred on X, T = 15', 6', about half-way between 
the 5 kpc and the 10 kpc arms did not show emission above 
0.13Kkms^' (1 s.d. at 23" resolution), which is 30 times be¬ 
low the average intensity along the corresponding segment at 
the crest of the 10 kpc arm. 

The apparent half-power width of the arms on Fig. Eb is 1-2 
kpc. This includes the true arm width within the plane of M31, 
the arm thickness perpendicular to the plane and the effect of 
beam-smoothing. Nieten J2001t attempted to disentangle these 
3 effects for the NW part of the 10 kpc arm and derived a ‘true’ 
width of 500 + 100 pc and an arm thickness of 150 + 50 pc. 

We have repeated the above analysis for the Hi map of 
Brinks & Shane (HUS, smoothed to the same resolution of 
45", for the 20 cm radio continuum map of Beck et al. (I1998t 
at the original 45" resolution and for the 175 jjm map at 90" 
resolution (see Fig. Et). The arm-interarm contrast of the in¬ 
tegrated Hi-line brightness is ^ 4, which is 5 times smaller 
than that of CO in the same disk section, and the apparent half¬ 
power arm width of 6' is 3 times larger. Contrary to CO emis¬ 
sion, H I emission is detected everywhere between the arms at 
5 kpc and 10 kpc. At 420 cm the contrast is about 2.5 and the 
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Radius [kpc] Radius [kpc] 


Fig. 4. (a) A family of logarithmic spiral segments (defined in the plane of M 31) fitting the CO arms in the NW quadrant, 
superimposed onto the velocity-integrated CO emission, /i_o. smoothed to a resolution of 45". The spiral segments densely 
cover the NW quadrant between Y — V and Y - 14' on the minor axis (indicated by the arrow), (b) Top: Mean profile of the 
CO intensity across the 5 kpc and 10 kpc spiral arms, averaged along the spiral segments shown in Fig.|3t. Middle and Bottom: 
same as the upper figure for the SW, SE and NE quadrants. The scale on the abscissa shows the galactocentric distance along 
the minor axis (where 1' corresponds to 1050 pc) for the NW and SE quadrants, and along the major axis for the SW and NE 
quadrants (where 1' corresponds to 227 pc). The arrows indicate the points of intersection of the two-armed logarithmic spiral 
pattern with the minor or major axis; note that they do not always coincide with the CO arms, (c) Comparison of the mean CO, 
HI, 20 cm continuum and 175 pm continuum emission profiles in the NW quadrant. The H i data of Brinks & Shane d 19841 1 were 
smoothed to an angular resolution of 45", the original resolution of the 20 cm continuum data of Beck et al. (I1998t . The 175 pm 
continuum data of Haas et al. dl998> have a resolution of 90". The CO profiles in the upper figure were obtained from the CO 
map smoothed to resolutions of 45" (full line) and 90" (dotted line), respectively. 


arm width is in between the widths in CO and H i. This was 
also noted by Berkhuijsen et al. (IT^ at a resolution of 75", 
who found that the width at A20 cm corresponds to the width 
of the total gas arm. 

The molecular arms traced by CO are much narrower and 
thinner than the H i arms and the arm-interarm contrast is much 
higher. This indicates that the molecular phase is short-lived 
compared to the life time of the H i gas. 

Erom a comparison of CO emission and visual extinction 
in a wide strip centred on the SW bright arm Neininger et al. 
GM) concluded that CO(l-O) line emission is a good tracer 
of the molecular gas in M31, including the interarm region. 
The > 5 times larger arm-interarm ratio in CO than in H i thus 
implies that the molecular gas has almost vanished in the inter¬ 
arm regions. Yet, molecular and atomic clouds have about the 
same velocity dispersion (see Sect. 5), have roughly the same 


response to the stellar potential and follow the same orbits with 
the same orbital velocity. This indicates that molecular clouds 
become mostly atomic when leaving the arms. 

An accurate determination of the molecular cloud lifetimes 
would require a density-wave model of M 31 that explains the 
observed CO arm pattern. Such a model is outside the scope 
of the present paper. A crude upper limit to this lifetime can 
nevertheless be estimated from the CO and H i arm-interarm 
ratios. Consider a gas cloud orbiting around M 3rs center and 
crossing a spiral arm. Erom Fig. Et we estimate that the total 
gas arm-interarm ratio is 4-5. This means that, due to stream¬ 
ing motions, the clouds stay 4-5 times longer inside the arms 
than they would if they followed purely circular orbits at a con¬ 
stant velocity. For an arm width of 0.5 kpc, a pitch angle of 8° 
and an orbital velocity of 300km s ', the time spent in an arm 
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is 5 10^ yr, which is an upper limit to the age of the molecular 
clouds. 

3.3. Line ratios in seiected regions 

Our observations of several selected regions on a hner grid 
yielded maps with high signal-to-noise ratio in the (2-1) line as 
well as in (1-0). Two example maps of the velocity-integrated 
CO (2-1) line intensity (Fig. 0 cover bright segments, about 
2-A kpc long, of spiral arms: one located near the major axis 
about 6 kpc from the centre (Fig. Bi), the other west of the 
main axis and about 9 kpc from the centre (Fig. 12)). The ratio 
f? 2 i of the (2-l)-to-(l-0) line intensities is often regarded as 
a hrst indicator of excitation conditions in the clouds emitting 
the CO lines. To apply this test to our data, we smoothed the 
(2-1) map to the same resolution (23") as the (1-0) map, se¬ 
lected those points were both line intensities were above 3cr, 
and computed the line ratio. In the region of Fig. Bi we hnd 
an average ratio of 0.65 + 0.1 with values reaching up to 1.0 
near D615 {X,Y = 25:5, 3;0);^ in the clouds of Fig. |5t we 
hnd an average of 0.5 ±0.1 with ratios up to 0.7 around D84 
(X,T = -16:5, -8:7). 

These values are compatible with simple standard assump¬ 
tions about the CO excitation, i.e., that (i) both lines are op¬ 
tically thick, (ii) they have the same excitation temperature 
Tex which equals the gas kinetic temperature and (iii) 
that they sample the same gas of uniform conditions. The ex¬ 
pected line ratio under these conditions is /? 2 i = 0.49 for 
r = Tex = Tkin = 3.5 K, R 21 = 0.67 for T = 6 K, and 
R 2 \ - 0.79 for T = 10 K. The line ratios we hnd in the spi¬ 
ral arms of M 31 are also similar to values found widely over 
the Milky Way (e.g., Sakamoto et al. ll997t as well as those of 
0.7 to 0.8 seen in M51 (Garcia-Burillo et al. ITWn and 0.60 
to 0.85 in NGC 891 (Garcia-Burillo et al. [T^ . For these rea¬ 
sons the CO excitation in the bright regions of our survey is 
clearly different from that of gas which shows very low “sub- 
thermal” values of f? 2 i ~ 0.2 to ~ 0.4 found, e.g., by Allen & 
Lequeux (I1993t and Loinard et al. J 19951 towards some posi¬ 
tions in M 31. Of course these papers concern mostly positions 
of small clouds or weak emission, whereas our values are se¬ 
lected for positions with relatively strong lines and therefore 
probably high ratios /? 2 i- We note, however, that Melchior et 
al. (I2000t reported an /? 2 i ratio of 0.65 towards the weak CO 
complex associated with D395A, located only 350 pc from the 
centre of M 31. This value is quite similar to the ratios that we 
hnd. A more complete study of CO excitation in M 31 would 
have to include positions of weaker emission and also observa¬ 
tions of optically thin CO isotopes; this is beyond the scope of 
our present paper. 


Y offset in kpc 

4 2 0 -2 -4 



Y offset in kpc 
-7 -8 -9 



Y offset in arcmin 

Fig. 5. Two maps showing the integrated emission of the 
'^CO(2-l) transition at a resolution of 12". The contour lines 
correspond to 3, 5, 7, .. .times the noise scale), (a) Spiral 
arm on the major axis at a distance of 5.5 kpc north of the cen¬ 
tre; the noise is 0.7 K km s ^ (b) Spiral arm region about 9 kpc 
west of the centre; the noise is 1.2 Kkm s“*. 


^ The dark cloud complexes (e.g. D615) are indicated with their ap¬ 
pellation in Hodge’s Atlas of the Andromeda Galaxy (Hodee ll981a> . 
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M31 CO HPBWz23'' North (superposed on U band) 



45- 30’ 15’ 00’ 


M31 CO HPBWz23" South (superposed on U band) 



-30’ 
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Fig. 6. Contours of '^CO( 1-0) emission from M 31 at 23" resolution on a U-plate of Walterbos & Kennicutt J1988t . The contour 
levels are 1 (= 3x r.m.s. noise), 4 and 8 Kkm s“^. (a) Northern half of M 31. (b) Southern half of M 31 
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4. Comparison of CO with other gas and dust 
tracers 

The high resolution data presented in this paper allow detailed 
comparisons with other tracers of the gas and the dust. Figure]^ 
shows contours of CO brightness superimposed onto a U-band 
image (Walterbos & Kennicutt [MHli. The CO closely traces 
the dust lanes, especially in the western part of M 31 (T < O') 
where the dust lanes stand out against the light of the stellar 
bulge. The extended dark region around X = 30', Y - -7' 
shows little CO, but this area has low brightness in many con¬ 
stituents of M 31 (see Fig.l^l. The weak inner CO arm crossing 
the major axis near X = -12' appears to coincide with a narrow 
dust lane seen at 7.7 fim by Willaime et al. J2001t who observed 
the central part of M 31 with ISOCAM LW6 (see Fig.0i. The 
CO clouds near the major axis between X - 0' and 10' also 
coincide with dust features at 7.7 fj.m, but not all dust features 
are seen in CO. 


M31 EOCAM LW6 + PICOVELETA CO (1-0) 



Fig. 7. Contours of ^^CO(l-O) emission from M 31 at 23" res¬ 
olution on the dust emission at Al.1 fim observed in the cen¬ 
tral bulge (Willaime et al. 1200 II Fig. 2). The contour levels are 
1 (= 3x r.m.s. noise), 4 and lOKkms”'. The white dot near 
the centre of the image indicates the centre of M 31. Note the 
alignment of CO clouds along the faint elliptical ring south of 
the nucleus. 

Figure |3 collects the observed brightness distributions of 
CO, H I, total neutral gas, FIR at 175 /tm and 20 cm radio con¬ 
tinuum. The 175 jum emission mainly traces cold dust at tem¬ 
peratures near 16 K (Haas et al. fT^ . The difference between 
the distributions of CO and Hi is striking. The Hi arms are 
much smoother and wider than the CO arms and weak Hi 
emission is seen nearly everywhere outside the spiral arms. 
We obtained the distribution of the total gas column density, 
Ngas = N(H i) + 2A(H2), after smoothing the CO map to the 


angular resolution of 24" x 36" of the Hi map (Brinks & 
Shane [1984). We used the Galactic conversion factor Aco = 
A^(H 2 )// r„b('^COi-o)dv = 1.9 X 1020molcm-2(Kkms-i)-' 
given by Strong & Mattox J1996> . which was assumed to be 
constant across the CO map. This conversion factor is sup¬ 
ported by 1.2 mm observations of the thermal emission from 
dust in two extended regions in the disk (Zylka & Guelin, 
in prep.) and by virial mass estimates of several molecu¬ 
lar cloud complexes in M31 (Muller Guelin et al. 

The Hi column density was calculated from the rela- 
tionA(Hi)/ jTbdv = 1.82 10^® at cm“^, valid for optically thin 
lines. This relation may lead to a significant underestimate of 
A(H i) if Tb > 70 K (Braun & Walterbos ll992t . i.e. at the crest 
of the arms. The distribution of the total gas calculated this way 
is shown in Fig. HJi. 

Comparing the various distributions in Fig. we notice 
that the nuclear area is only prominent at 175 /tm and A20 cm 
radio continuum, but all distributions in Fig. |3 show the pro¬ 
nounced ring of emission at about 10 kpc from the centre, 
where also most of the star-formation regions are located. The 
spiral arm closer to the nucleus, clearly visible in CO and at 
175jum, is hardly seen in Hi and radio continuum. This in¬ 
dicates that the CO-to-Hi brightness ratio increases towards 
smaller radii. Only the Hi and the 175/im distributions show 
extended weak emission in between the spiral arms. Thus the 
CO emission is concentrated to regions of the denser clouds 
seen at 175/im and in Hi, located in the spiral arms, and does 
not trace weak and extended interarm emission at 175/im, 
which is especially visible at T > O'. In Sect. 4.2 we present 
first results of a correlation study between CO, H i and 175 /im. 
Radial distributions and the gas-to-dust ratios are shown in 
Sect. 4.1. 

The radio continuum emission at A20 cm (Beck et al. ll998> 
is also concentrated in the main spiral arms that form the emis¬ 
sion ring. As at this wavelength most of the emission is non- 
thermal, magnetic fields and cosmic rays are concentrated in 
this ring. Berkhuijsen ( 119971 . using the CO data of Koper et al. 
itTWni . compared nonthermal emission with CO and H i data 
at a resolution of 9'. Although significant correlations indeed 
exist at this resolution, detailed correlations with the new data 
are required to enable an interpretation in terms of the coupling 
between magnetic fields and gas. This will be the subject of a 
forthcoming study. 

4.1. Radial distributions 

In Fig.|9lwe show the radial variations of the observed CO and 
Hi brightness distributions averaged in 230 pc-wide circular 
rings in the plane of M 31 (i.e. 1' on the major axis). Because 
the spiral arms in M 31 are tightly wound, the narrow CO arms 
clearly stand out in the radial distributions. In the northern part 
(X > O'), the average intensity of the inner arm at R ^ 25' is 
comparable to that of the arms forming the bright emission ring 
at R ^ 50', whereas in the southern part (X < O') this inner arm 
is much weaker. The radial profiles of H i are much wider than 
those of CO, but the inner arm and main ring are still indicated. 
The HI disk extends as far as R ^ 150', 34 kpc from the nucleus 
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(Emerson 1 1974> . about 2 times further than the molecular gas 
in the map of Koper et al. (I1991> . 
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Fig. 8. Distributions of neutral gas, cold dust and radio continuum in M 31. From top to bottom: (a) emission observed in the 
'^CO(l-O) line (this paper), (b) emission observed in the Hi line (Brinks & ShaneQ^HD, (c) emission from the total neutral gas, 
N(Hi)+ 2N(H2), with Xco = 1-9 x 10^*^molcm“^(Kkms“')"' (see text), (d) emission from cold dust at AllS/im (Haas et al. 
119981 . (e) radio continuum emission at A20 cm (Beck et al. ll998t . The half-power beamwidth is indicated in the upper right-hand 
corner of each map. The intensities at the maximum near X,Y - -42', 8' (white-red) and the minimum near X,Y - -30', 10' 
(blue-green) in each map are: (a) 6.7 and 0.8 in 10^' atcm^^, (b) 6.5 and 2.0 in 10^^ atcm^^, (c) 12.1 and 2.6 in 10^' atcm^^, 
(d) 85 and 27 in MJy/sr, (e) 2.9 and 0.7 in Jy/beam area. The M 31 coordinates are based on the centre position (a, (5)so - 
(0H0™0(k3, 41°00'03") and a position angle of the major axis PA - J>T.l (see footnote 1). 
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Our CO profiles of the southern half are only half as strong 
as those shown in Fig. lOd of Loinard et al. J1999t . whereas the 
HI profiles are nearly identical. This may indicate a difference 
in the adopted CO intensity scales, e.g. in the calibration or in 
the corrections for beam efficiencies. We note that in Sect. 4.3 
we find a close agreement of our CO intensity scale with that 
of the 1.2-m telescope at the Center for Astrophysics (Dame 
et al. I1993> . to which the standard Galactic Xco is calibrated. 
We also find that even after smoothing our map to the effective 
resolution of 1' of the map of Loinard et al., the half-intensity 
width of the CO arms in their Figs. 5 and 6a is about 1.5 times 
larger than in our map. We may attribute this discrepancy to the 
much higher error beam of the FCRAO telescope compared to 
that of the IRAM 30-m dish. 


0 


5 


radius [kpc] 

10 15 


20 



Fig. 9. Radial profiles of /i_o (full lines, left-hand scale) and 
/hi (dashed lines, right-hand scale) for the northern half (N), 
the southern half (S) and the full area of M 31 (Nh-S). The pro¬ 
files show intensities along the line of sight averaged in circular 
rings of 1' width in the plane of M 31; typical standard devia¬ 
tions are 0.01 Kkms ' and 5 Kkms ' for the CO and the Hi 
profiles, respectively. The CO profiles were obtained from the 
map at 23" resolution (Fig. [3 and the H i profiles from the map 
of B&S at 24" X 36" resolution. The limited extent of the CO 
map on the major axis, R - 55' in the north and R - 50' in 
the south, reduces the intensities of the profiles of /i_o beyond 
these radii by a few percent. 


The fraction of molecular-to-total gas, 2A1(H2)/Algas, gen¬ 
erally decreases with increasing distance from the centre (see 


M31 2N(Hp/N,„ 



Fig. 10. Radial variation of the molecular gas fraction 
2A1(H2)/Algas for the northern half (N, upper left-hand scale), 
the southern half (S, right-hand scale) and the full area of M 31 
(N-hS, lower left-hand scale), where Agas = A(Hi) + 2A(H2) 
and Xco - 1.9 x 10^°mol cm“^ (Kkms^')^' (Strong & 
Mattox [T996> . The CO map was smoothed to the resolution of 
24" X 36" of the H i map before the molecular fractions were 
computed. The molecular fractions were calculated from the 
column densities averaged in circular rings of 1' width in the 
plane of M 31. 


Figs. 1101 and llH . The highest values occur on the weak inner 
arm at /? ^ 12' (2-3 kpc) and on the bright inner arm atR ^ 25' 
(5-6 kpc) where the molecular fraction is nearly 0.5 in the north 
and about 0.25 in the south (see Fig. Go). On the bright ring at 
R ^ 50' the molecular fraction is about 0.2 both in the north¬ 
ern and in the southern part, although the northern part of the 
ring contains about 1.5 times more neutral gas than the south¬ 
ern one. The averages in circular rings underestimate the ratios 
on the arms. Figs. Et and yield a decrease from ~ 0.6 on 
the arms at 5 kpc to 0.3-0.4 on the arms at 10 kpc. The ob¬ 
served values of 2 A(H 2 )/Agas and its radial decrease are typi¬ 
cal for nearby galaxies (Honma et al. I1995t . The decrease of 
the molecular fraction in Fig. [lOl confirms the decrease in the 
fraction of molecular gas mass along the major axis of M 31 
reported by Dame et al. (I1993> . Along the arms the molecu¬ 
lar fraction varies considerably (see Fig. 11 1 > . The highest frac- 
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M31 2N(H2)/Nga3 HPBW=24" x 36" 
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Fig. 11. Distribution of the molecular gas fraction 2A(H2)/A^gas in M 31 at a resolution of 24" x 36". Only data points larger than 
3x r.m.s. noise in the CO and H i maps were used. Agas = A^(H i) + 2A(H2) and a constant value of Xco = 1.9 x 10^° mol cm^^ 
(Kkms''')“* (Strong & Mattox [T996> were assumed. 


tion detected is 0.96 for a cloud near the northern major axis at 
X — 13', Y = -1', which seems purely molecular. 

We smoothed the distributions of CO, Hi, total gas and 
20 cm radio continuum in Fig. |3 to the angular resolution of 
90" of the 175 fim map and compare their radial distributions 
in Fig.[21 In all constituents the main ring atR ^ 50' is brighter 
in the north than in the south. The pronounced molecular arm 
in the north at ^ 25' is invisible in radio continuum, possibly 
because of a lack of relativistic electrons (Moss et al. 119981 . 
Disregarding the central region R < 10', the profiles of molec¬ 
ular gas and 175 /im emission from cold dust are most alike. 

Using the profiles in Fig. El we calculated the apparent 
gas-to-dust ratios 2 A(H 2 )//i 75 , A^(Hi)//i 75 and A^ga.s//i 75 as a 
function of radius which are presented in Fig. El The ratio 
2 A(H 2 )//i 75 is clearly enhanced in the spiral arms at ^ 25' 
and R ^ 50', especially in the north, whereas A(H con¬ 
tinuously increases steeply from the centre outwards by nearly 
a factor 20. As H i is the dominant gas component, the appar¬ 
ent total gas-to-dust ratio A^gas/fns increases by about a factor 
20 between R - 0’ and R - 60'. However, the physical real¬ 
ity of this strong increase may be questionable. Firstly, the Hi 
line opacity changes with radius; secondly the Tfco conversion 
factor may vary with radius; thirdly, Im may not reflect the 
dust column density if the dust absorption cross section and, 
especially, the dust temperature vary with radius. 

The first two causes are unlikely. Braun & Walterbos J 19921 
showed that the atomic gas temperature in M 31 increases to¬ 
wards larger radii, while the 21-cm line opacity decreases. This 
variation, however, is less than ~ 20 % and cannot account for 
the strong gradient in A(H i)//i 75 . 

The behaviour of Xco with radius has been the subject of 
several studies in the Milky Way and nearby galaxies (Wilson 
119951 Sodroski et al. 119951 Strong et al. 120041 . Xco is found 
to increase with increasing radius, perhaps in relation with de¬ 


creasing metallicity. Such an increase, if present at all in M 31, 
would only enhance the radial variations in A^(H 2 )//i 75 and 
iVgas/7l75- 

A study of 50 bright cloud complexes at a resolution of 
at least 12 ", which is high enough to alleviate rotation veloc¬ 
ity gradients, at radial distances between 5 and 12 kpc (Muller 
l2^ Guelin et al. lMHll yield Xco = (1 - 5) lO^^mol cm^^ 
(Kkms ') * without a radial dependence. Melchior et al. 
J20001 and Israel et al. <19981 find similar values for two com¬ 
plexes at R = 0.35 kpc and R = 2.5 kpc. Hence, except per¬ 
haps for the observation by Allen & Lequeux <19931 of a dark 
cloud (D268, R = 2.5 kpc) with low CO luminosity and large 
linewidth, which may not be in equilibrium, there is no evi¬ 
dence of a radial variation of Xco in M 31. 

It is more difficult to rule out a variation of the dust tem¬ 
perature Ti with radius. The thermal emission of cold dust at 
175 jum varies like Td to the power 4-5, so even a mild decrease 
of Td outwards could cause a strong decrease of / 175 . In fact, 
the apparent increase of Xgas//i 75 from 1 to 14 kpc could be 
explained by a decrease of the dust temperature from ~ 25 K to 
15 K. Only extensive mapping of this emission at longer wave¬ 
lengths will allow to discriminate between a temperature effect 
and a dust-to-gas variation. 

An increase of the apparent gas-to-dust ratio with increas¬ 
ing radius in M 31 was found by several authors from com¬ 
parisons of A(Hi) and optical or UV extinction (Bajaja & 
Gergelv lT977l Walterbos & KennicuttESl Xu & HelouESl 
Nedialkov et al. 120001 Savcheva & Tassev 120021 . The latter 
authors used globular clusters and found a similarly strong in¬ 
crease of the apparent gas-to-dust ratio as we do. The distri¬ 
bution across M 31 of the reversed ratio, i.e. the apparent dust- 
to-total gas ratio / 175 /Xgas, is shown in Fig. El The general 
decrease from the centre outwards is clearly visible. / 175 /Xgas 
tends to be low on the gas arms where the CO emission is 
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M31 50’N 0’ 50’S 



Fig. 13. Radial variation of apparent gas-to-dust ratios for 
the northern (left panel) and the southern half (right panel) 
of M31 derived from the profiles in Fig. fT^ Thin full line 
- Ngas/Ins, dashed line - A(Hi)//i 75 and thick full line - 
2A(H2 )//i 75. A^gas = A^(Hi) + 2A(H2) and a constant value of 
Xco = 1-9 X 10^° mol cm“^ (Kkms^')^' (Strong & Mattox 
[T9^ were used. Note the increase of 2 A(H 2 )//i 75 on the 
molecular arm at R ^ 25' and on the main emission ring at 


Fig. 12. Radial variations for the northern (left panel) and the 
southern (right panel) half of M 31 of A^gas, A^(H i) and 2 A(H 2 ) 
(left-hand scale) as well as of Insum (inner right-hand scale) 
and 7(20 cm) (outer right-hand scale) obtained from the distri¬ 
butions in Fig. 1^ The gas maps and the radio continuum map 
were smoothed to the resolution of 90" of the 175 ;um map. 
The line-of-sight surface brighnesses were averaged in circular 
rings around the centre of 5' width in the plane of M31. The 
two arrows indicate the extent of the CO map on the major axis. 

strong. Even the many roundish minima in the distribution of 
/ 175 /Agas coincide with CO clouds. 

It is interesting to note that on the brightest spot of emission 
at 175 ;um, about 4' west of the nucleus, the value of / 175 /A^gas 
is normal for that radius. This bright emission comes from one 
of the coldest dust clouds in M 31 analyzed by Schmidtobreick 
et al. (I2000t . Willaime et al. ll2()()l> showed that this cloud is 
located at a junction or superposition of several thin dust lanes 
seen at 7.7 //m (see Fig.Q, which may explain its normal dust- 
to-gas ratio. 

4.2. Correlations between the distributions of CO, H 1 
and FIR 

In Sect. 3.2 we compared the CO and Hi distributions across 
the spiral arms. We now compare the general distribution of 
CO with that of Hi and each of them with the 175/im FIR 
emission using the CO and H i maps smoothed to 90", the res¬ 
olution of the FIR map. At this resolution the CO and H i arms 
near R - lOkpc merge to a broad emission ring, similar to that 
seen at 175 fim. The comparison was restricted to radii R < 60' 
(= 13.6 kpc) and intensities above 3x the r.m.s. noise at 90" 


R 50'. 


resolution. Therefore, the noise in the CO map determined the 
selection of points in correlations involving CO intensities and 
the noise in the H i map that in correlations involving intensities 
of Hi or total gas. In spite of this difference, the selected pix¬ 
els are largely the same; most pixels are located on the broad 
emission ring. To obtain sets of independent data points, i.e. 
a beam area overlap of < 5%, only pixels spaced by > 1.67x 
the beamwidth were considered. The resulting correlations be¬ 
tween CO, Hi and the FIR are listed in Table Inland examples 
of correlation plots are shown in Fig. 1151 

In view of the distinct differences between the CO and H i 
distributions, it is not surprising that their velocity-integrated 
intensities are not well correlated (see Fig. [T3t and b). In the 
panel showing the full radial range, 7? = 0 ' - 60', two maxima 
occur which correspond to the inner bright CO arm and the 
main emission ring, respectively. In the inner arm the molec¬ 
ular fraction is much larger than in the main ring, and the de¬ 
crease of this fraction outwards obviously contributes to the 
large spread in the points. Therefore we also analyzed the inner 
(R - 0' -35') and outer (R - 35' -60') radial ranges separately. 
In the inner part 7i_() and 7 hi are not correlated (the correlation 
coefficient p - 0.3, see Table |2j, but for the outer range some 
correlation is visible (Fig. 115b). with p - 0.61 +0.09.Note that 
the correlation follows a power law with exponent 1.7 + 1.0. An 
exponent of 2 is expected if the formation and destruction rates 
of H 2 are balanced (Reach & Boulanger ll998t . This occurs in 
the denser H i phase at temperatures near 80 K, which are also 
found in M31 (Braun & Walterbos ll992> . For a better check 
of this possibility in M 31 this exponent could be determined in 
narrow radial ranges at the original resolution of 24" x 36" of 
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Fig. 14. Distribution of the apparent dust-to-total gas ratio, / 175 /A^gas (MJy sr ')/(10^^ at cm ^), in M 31 at 90" resolution. Local 
minima (holes) are caused by enhancements in A^(H 2 ) (compare Fig. m or Fig. CHi. 

Table 2. Correlations between CO, H i and 175 /tm 


X 

Y 

AR 

(') 

log Y = a log A' + 
a b 

Correl. 
coeff. p 

N" 

/hi 

h-Q 

0-35 

2±5 

-5+20 

0.3+0.1 

90 

(Kkms^') 

(Kkms^*) 

35-60 

1.7+1.0 

-5.2±2.4 

0.61+0.04 

234 



0-60 

1.3+0.5 

-4.1 + 1.4 

0.50+0.04 

324 

/l75 

^1-0 

0-35 

1.6+1.1 

-2.53+0.63 

0.56+0.07 

94 

(MJy sr“') 

(Kkms^*) 

35-60 

1.57+0.20 

-2.39+0.16 

0.81+0.02 

234 



0-60 

1.60+0.22 

-2.47+0.16 

0.75+0.02 

328 

/l75 

^Hl 

0-35 

0.8+1.8 

1.4+2.2 

0.51+0.07 

125 

(MJy sr-‘) 

(Kkms^*) 

35-60 

0.84+0.27 

1.85+0.44 

0.75+0.03 

264 



0-60 

1.09+0.32 

1.36+0.42 

0.54+0.04 

388 

/l75 

Afgas 

0-35 

1.42+0.53 

0.89+0.60 

0.68+0.05 

134 

(MJysr-‘) 

(10'* atcm^^) 

35-60 

1.12+0.14 

1.64+0.22 

0.85+0.02 

263 



0-60 

1.23+0.19 

1.27+0.27 

0.67+0.03 

395 


'^Weighted fits of the bisector through N pairs of (log A, log Y), where N is the number 
of independent points (Isobe et al. 1 19901 Nieten l^OH 


the HI map, which may yield better correlations than obtained 
here. 

The middle panels of Fig. [21 show that in the interval R = 
35' - 60' both /hi and /i_o are well correlated with the dust 
emission at 175 /tm. With correlation coefficients of p = 0.75 + 
0.03 and p - 0.81 + 0.02, respectively, both correlations are 
highly significant (Table|2j. The relationship between /hi and 
/i 75 is about linear, whereas that between /i_o and /175 follows a 
power law with exponent 1 . 6 + 0.2 (see Table|2j- This difference 
reflects the power-law dependence of /i_o on /hi in Fig. [BJ). 
The two branches visible at log/i 75 < 1.3 in Fig. ini; are caused 
by the strong radial gradient in the apparent atomic gas-to- 
dust ratio A^(Hi)//i 75 plotted in Fig. Il3] points in the upper 
branch are at larger radius than points in the lower one. This 


branching does not occur in the correlation between /i_o and 
hi5 (Fig. I15I1'). because 2N{\i2)lhi$ increases much less with 
R than A(H i)lI\i$. 

Correlations between the total gas column density and the 
175 pm emission for both the inner and the outer part of M 31 
are better than those between each of the gas components alone 
and /i 75 (see Table|2land Figs. [Bt and f). The correlation for 
the outer part is very good (p = 0.85 ± 0.02), much better than 
for the inner part (p = 0.68 + 0.05). Because in the inner part 
A^gas is dominated by the molecular gas, but in the outer part 
by the atomic gas, the power-law exponent in the inner part 
(1.42 + 0.53) is larger than in the outer part (1.12 + 0.14). The 
latter exponent is consistent with a linear correlation between 
A^gas and / 175 . Linear correlations between gas and dust in M 31 
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were also found by other authors. Pagani et al. ( fT^ found 
a linear dependence of the intensity in the ISOCAM LW2 fil¬ 
ter (5.0-8.0;um) on Ng^s in ^ field centred on the southern ma¬ 
jor axis near X - -30'. In a spiral-arm region centred near 
X - -20', Y = -7' in the SW quadrant, Neiningeret al. raj 
obtained a good linear correlation between Agas and the ap¬ 
parent red opacity. In the Milky Way gas and dust were found 
to correlate in detail (Bohlin et al. 119781 Boulanger & Perault 
fT988l Boulanger et al. lT^ . 

One may ask why the 175 fim dust emission follows more 
closely Agas than A(H i) alone, since the gas is mainly atomic 
(see Table|3 see also Buat et al. ll989t . There are several rea¬ 
sons for this. First, emission from dust in Hi clouds as well 
as from dust in H 2 clouds contributes to the 175/im emis¬ 
sion. Second, according to dust models (see e.g. Mezger et al. 
[I990j, the dust absorption coefficient and the dust emissivity 
per H-atom are two (or more) times larger in the dense molec¬ 
ular clouds than in the diffuse atomic clouds. Hence, for equal 
dust column densities and temperatures, the emission from dust 
in CO clouds exceeds that of dust in Hi clouds as soon as 
2A(H2)/A(Hi) > 0.3. 

The data presented in Sects. 4.1 and 4.2 show that a ra¬ 
dial gradient in the apparent gas-to-dust ratio exists that in¬ 
fluences the correlations between the gas components and the 
dust. Corrections for H i-line opacity, better knowledge of the 
conversion factor Xqo and, probably most important, a deter¬ 
mination of the dust temperature as a function of radius are 
needed to understand whether this radial gradient reflects a real 
increase in the dust-to-gas ratio in M 31. In the meantime com¬ 
parisons of gas and dust on small scales in narrow radial ranges 
would be of interest. 

4.3. Molecular and total gas mass 

The CO intensity integrated over the area of M 31 gives an es¬ 
timate of the total molecular mass if one assumes a conversion 
factor between Ii-o and A(H 2 ). We integrated the CO inten¬ 
sities in Fig. out to a radius of 80' in the plane of M31, 
which corresponds to R - 18 kpc. With the same conver¬ 
sion factor Aco as used above we And a molecular mass of 
^(Ha) = 3.45 X 10* Mq. Along the major axis our map only 
extends to R 55', so we missed the emission near the major 
axis between R = 55' and 75' detected by Koper et al. ( fTWn . 
From Figs. 3b and 10c in Dame et al. ( fT^ we estimate that 
this contributes about 5% of the total. After correcting for these 
5%, we get M(H 2 ) = 3.6 x 10* Mq within a radius R - 18 kpc. 
This is in excellent agreement with the value obtained by Dame 
et al. (I1993> . which is M(H 2 ) = 3.45 x 10* Mq after scaling it 
to the distance of 780 kpc adopted here. Given that Dame et al. 
used the same value for Aco as we do, we can also conclude 
that their CO intensity scale of radiation temperature, Tr, is in 
close agreement with our scale of main beam brightness tem¬ 
perature, Tnib. 

To obtain the atomic gas mass within R = 18 kpc we inte¬ 
grated the Westerbork H i map of Brinks & Shane J1984> giv¬ 
ing M(Hi) = 2.6 X 10^ M 0 in the optically thin case. Thus 
for R < 18 kpc the H 2 mass is 14% of the Hi mass and 12% 


of the neutral gas mass. The total H i mass of M 31 at the dis¬ 
tance of 780 kpc is 4.86 X 10^ Mq (Cram et al. ll980l corrected 
by Dame et al. I1993t which gives a total neutral gas mass in 
M31 of Mgas = 5.2 X 10^ Mq and a molecular mass fraction 
M(H 2 )/Mgas = 0.07. 

The bulk of the dust in M31 is cold dust at a tempera¬ 
ture of Td = 16 K. Haas et al. J1998t estimated a total dust 
mass of Md = 3.8 x 10^ Mq at D - 780 kpc. However, 
Schmidtobreick et al. J2000> found Md = 1.3 x 10^ Mq using 
the same data but a different method of calculation. Hence the 
apparent gas-to-dust mass ratio in M31 is Mgas/Md = 137- 
410. If the HI masses are corrected for opacity of the H i, they 
increase by 19% (Braun & Walterbos ll992t and the ratio be¬ 
comes Mgas/Md = 163—488. For the Milky Way Sodroski et al. 
raj obtained a gas-to-dust mass ratio for the entire Galaxy 
of Mgas /Md - 160 + 60, similar to the lower estimates for M 31. 

5. The CO velocity field 

Figure shows the first moment of the CO spectral data cube, 
i.e. the intensity-weighted mean CO velocity as a function of 
position. This figure can be compared with Fig. 15a of Brinks & 
Burton ( 119861 - hereafter B&B), which shows the first moment 
of the Hi data cube. Contrary to Hi, CO emission becomes 
very weak at galactocentric radii R > 18 kpc, where the Hi 
disk starts to flare and to warp (scaled to the distance of M 31 
of 0.78 Mpc). The velocity held of Fig.^ is therefore little af¬ 
fected by foreground and background emission from the warp 
and is more appropriate to derive the velocity held in the disk 
than the velocity held of H i. Moreover, as M 31 is at a galac¬ 
tic latitude of b - -22°, contamination by foreground Milky 
Way clouds at velocities of interest here (< -50 km s“') is un¬ 
likely. The average deviation between the Hi and CO mean 
velocities is only ^ 10 km s ', which may imply that the bulk 
of the HI emission in the region covered by our CO map arises 
from M31’s disk, although the different velocity resolutions 
(8.2km s ' for Hi and 2.6kms“' for CO) may contribute to 
the differences. 

The distribution of the mean CO velocities is dominated 
by the rotation around the dynamical centre of the galaxy 
that has a velocity Vsys = -315kms“'. Only one small 
cloud near the centre of the galaxy and a large cloud com¬ 
plex near the northwestern minor axis show strong deviations 
from circular rotation. The small cloud is located at A, T = 
-OH, 3H (-250kms“'). Its non-circular orbit may be caused 
by the central bar tBerman r2001l Berman & Loinard l2002t . 

The cloud complex near the minor axis (A, Y = 3', -4'), 
that appears particularly bright at 175 fim, shows complex line 
profiles that may be related to the complex filamentary struc¬ 
ture visible in Fig. 0 Alternatively, they may trace streaming 
motions (Muller et al., in prep.). 

The width of the observed CO-line profiles varies strongly 
from position to position. The average linewidth (second mo¬ 
ment of the CO data cube) computed over the whole disk is 
10 + 5 km s“'; it is higher in the arms (15 + 5 kms^') and lower 
at the arm edges and in the interarm regions (8 + 4km s“'). 
The average CO linewidth is in agreement with the average 
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linewidth of the Hi emission of 8.1 km s' (at 1' resolution) 
arising from the disk (Unwin Il983> . 

Some of the line profiles are very broad (40-65 km s ', 
e.g. Fig. 12;) and show multiple velocity components. This is 
particularly the case for the cloud complexes associated with 
the dark clouds D47 and D39 (X,Y = -22(48, -7(53 and 
-41(9, 8(54, see Fig. |2^ and i) which show complex profiles 
with total widths of 50 and 65 km s ', respectively. We note 
that in the case of D 47 and D 39 the multiple profiles arise in 
the vicinity of bright H ii regions, but this possible connection 
is not a general feature. 

Other line profiles appear to be very narrow (< 5kms“'; 
see Fig. |2j>). The profiles associated with D 153 (X, Y - 
-17(8, -3(85, seeFig.Et) have a particularly high peak inten¬ 
sity combined with a narrow line width of only about 4 km s"'. 

A detailed analysis of the cloud-to-cloud velocity disper¬ 
sions inside the molecular cloud complexes, and of the inter¬ 
play between density-wave driven motions and peculiar mo¬ 
tions, linked to H ii regions or supernova remnants, will be the 
subject of a forthcoming paper (Muller et al., in prep.). 

The position-velocity diagram of the CO (1-0) intensity in¬ 
tegrated along the minor axis is shown in Fig. M The dominant 
features directly reflect the general velocity held in the disk 
(Fig.nh) and the fact that most CO emission is concentrated 
in spiral arms (Figs. □ 1 , 131 ). Due to the rotation of the M31 
disk, CO emission from the arms marked in Fig. Et is seen in 
a broad band of about 100 km s * width, stretching from X ^ 

-48', V -550kms"' to A" ^ h-55', v = -50kms“'. Most 

of the emission in this band is from regions near R = lOkpc. 
This large-scale picture agrees of course with that in the low- 
resolution survey by Dame et al. ifTWn . who also And some 
emission near X = -63, v = -550 km s ', outside the area 
of our survey. The CO spiral arms appear as loops within the 
range of this band (compare Braun J1991t for H i); most clearly 
seen in Fig. is the loop at the most negative X correspond¬ 
ing to the “dashed” spiral arm in Fig. Et Emission from the 
inner part of the “solid” spiral arm in Fig. Et is clearly visible 
at A" < 25', V ^ -200 km s“' above the main band of emission. 
Extended and relatively strong emission at velocities around 
-80 km s ' and X about 8' to 12' is from clouds located inside 
the “solid” spiral arm in Eig. Et, i.e. closer to the major axis 
(compare Eig.^i). While we cannot reliably trace spiral arms 
from the CO emission in this inner region, we note that these 
CO features agree in position and velocity with the inner loops 
of the spiral arms identified by Braun (1991) from Hi surveys. 
In a similar fashion, most smaller and weaker features at veloc¬ 
ities different from those of the main CO spiral arms in Eigs. Et 
an dEt are on other loops of Braun’s H i spiral arms. 

6. Summary 

The new '^CO(J=l-0)-line survey of the Andromeda Galaxy 
presented here covers an area of about 2° x 0(5, which is fully 
sampled with a velocity resolution of 2.6 km s ' and an angular 
resolution of 23", the highest angular resolution to date of a 
map of this extent. At the adopted distance of 780kpc the linear 
resolution is 85 pc x 400 pc in the plane of M 31. The On-the- 
Fly method of observing made it possible to measure nearly 1.7 


million spectra (before gridding) in about 500 hours of effective 
observing time on M31. The r.m.s. noise in the CO(l-O) line 
per 1 MHz channel is 25 mK in the northern half of the map 
and 33 mK in the southern half 

The velocity-integrated distribution, /i_o, and the velocity 
held are shown in Fig. □ The molecular gas is concentrated in 
narrow, filamentary arms between 4 and 12 kpc from the centre 
with a maximum near 10 kpc. The inner arm at ^ 5 kpc is 
remarkably bright. Only few clouds, often forming bridge-like 
structures, are detected in between the arms above 3x r.m.s. 
noise of typically 0.35 Kkm s '. The region within 2 kpc from 
the centre is almost free of molecular clouds brighter than the 
sensitivity of this survey. 

The thin CO arms define a two-armed spiral pattern that 
can be well described by two logarithmic spirals with constant 
pitch angle of 7° - 8°. At a resolution of 45" the arm-interarm 
contrast reaches a maximum of 20 in /i_o compared to 4 in 
Hi for the western bright arms. The Hi arms are much wider 
than the molecular arms, and diffuse H i exists everywhere in 
between the arms and at large radii. Few molecular clouds are 
visible outside a radius of 16 kpc. 

The velocity held of the molecular gas is very similar to 
that of the disk component in H i. At some positions perturbed 
velocity profiles occur that are possibly caused by nearby H ii 
regions. 

Several selected regions were also observed in the '^CO(2- 
1) line. At a resolution of 23" the line ratios are nearly constant 
with mean values of / 2 - 1 // 1-0 = 0.5 - 0.7 in the arms. These 
line ratios are similar to those observed in other galaxies and 
show no indication of subthermal excitation. 

Averaged radial profiles of the velocity-integrated CO and 
HI distributions show that for a constant conversion factor of 
Xco - 1.9xl0^''mol cm“^ (Kkms“')"'themolecularfraction 
of the neutral gas is enhanced on the spiral arms and decreases 
radially from about 0.6 on the inner arms to about 0.3 on the 
arms at ^ 10 kpc (see Fig. 11 H . Along the arms the molecular 
fraction also varies considerably. Comparisons with averaged 
radial profiles of the /1175 /im emission, which traces the cold 
(16 K) dust, revealed a strong, continuous increase of the ap¬ 
parent atomic gas-to-dust ratio of nearly a factor 20 between 
the centre and ^ 14 kpc, whereas the apparent molecular 
gas-to-dust ratio increases by about a factor of 4. The appar¬ 
ent total gas-to-dust ratio also increases by about a factor of 
20. 

The strong apparent gradients in the molecular fraction 
and the gas-to-dust ratios influence the cross-correlations be¬ 
tween CO, Hi and FIR(175 fim) intensities. In the radial range 
R - 35' - 60' (about 8-14kpc) the best correlation exists be¬ 
tween total neutral gas and 175 /im, followed by that between 
CO and 175 /urn and between H i and 175 fim. The relationships 
between Hi and 175 fim and total gas and 175pm are close to 
linear, whereas that between CO and 175 pm is a power law 
with exponent 1.6 due to a possible non-linearity between CO 
and H i. In the inner part of M 31, /? < 35' (but outside the nu¬ 
clear area), only total gas and FIR(175 pm) are reasonably well 
correlated (see TableEll. 

The total molecular mass of M 31 within a radius of 18 kpc 
is 3.6 X 10^ Mq, using the above-mentioned value of Tfco. As 
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Fig. 16. Position-velocity diagram along the major axis obtained by adding the intensities along lines parallel to the F-axis as 
function of position on the 2f-axis and velocity. The grey scale represents the intensities starting at contour level 0.5 K. The dotted 
line indicates the system velocity of -315 km s“'. 


the total Hi mass (without correction for opacity) is 4.86 x 
10® Mq, the total mass of the neutral gas is 5.2 x 10® Mq. The 
total mass of the cold dust is (1.3 - 3.8) x 10® Mq, hence the 
total gas-to-dust mass ratio is 410-137. The lower value is in 
agreement with the Galactic one. 

The wealth of information contained in our new '^CO(l-O) 
survey of M 31 enables a number of new investigations into the 
physical relationships between molecular gas, atomic hydrogen 
gas, cold and warm dust, ionized gas, relativistic electrons and 
magnetic fields. Such studies will be the subject of forthcoming 
papers. 

The data shown in Fig. [Qcan be obtained from M. Guelin 
(guelin@iram.fr). 
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